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ABSTRACT 

A detailed analysis was conducted for the AIMP-D spacecraft mis- 
sion to determine the launch conditions and trajectory shaping which 
would maximize the probability of attaining a lunar orbit with a lifetime 
in excess of six months. This study arose from the relatively large 
transfer- orbit injection errors which are associated with the spin- 
stabilized solid propellant injection motor. No midcourse correction is 
available; the only control parameter is the time of firing of the fixed- 
impulse lunar orbit injection motor. The large errors resulted in anal- 
ysis and solutions peculiar to this mission trajectory. It was found that 
for the applicable vehicle errors, maximum lunar orbit probability was 
obtained for a 72-hour flight time and a high injection flight path angle. 
A circular parking orbit could not be used because of the flight path 
angle requirement and the available launch times were thereby re- 
stricted. Analysis of the spacecraft and guidance constraints resulted 
in the selection of two launch periods, each of several days duration, 
occurring twice a year. 
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LUNAR ORBIT MISSION ANALYSIS FOR 
THE IMPROVED DELTA LAUNCH VEHICLE 
AND AIMP-D SPACECRAFT* 

by 

D. L. Mootchnik 
Douglas Aircraft Company, Inc . 
and 

J. Kork 

Goddard Space Flight Center 


INTRODUCTION 

Certain results of the Improved Delta AIMP-D Spacecraft Lunar Orbit Mission Analysis (Refer- 
ence 1) are discussed in this paper. The mission involves the placement of a payload consisting of 
an instrument package and a fourth stage retro motor into a translunar orbit by means of the three- 
stage Delta boost vehicle. In the vicinity of the moon, the spin- stabilized, fixed-impulse retro 
motor reduces the approach velocity, thereby causing capture of the spacecraft in a closed lunar 
orbit. 

Because of the high reliability of the Delta vehicle and the limited spacecraft weight capability, 
a basic management decision was made to minimize vehicle and guidance system hardware modi- 
fications. The addition of a midcourse correction system was precluded by this decision. Conse- 
quently, the trajectory errors introduced by the boost vehicle, particularly by the spin- stabilized 
solid propellant third stage, will propagate through the transfer trajectory and cause large, non- 
linear dispersions in the moon's vicinity. These potential large dispersions and the availability of 
only one control parameter, viz., fourth stage ignition time, have led to analysis methods and 
results unique to this lunar mission. 

The primary purpose of the mission analysis was to determine the launch conditions and 
trajectory shaping which would optimize the probability, for the given spacecraft weight, of obtain- 
ing a lunar orbit with a lifetime in excess of six months. Orbit lifetime is the period from lunar 
orbit insertion to the time when the spacecraft either escapes the moon's vicinity or impacts the 
lunar surface. Because of the trajectory dispersions and large orbit perturbations caused by the 


•This paper is based on work accomplished by the Douglas Aircraft Company, Missile and Space Systems Division, under Contract 
NAS7-265 (CCN 5) for Goddard Space Flight Center. 
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earth, a short lunar orbit lifetime is of primary concern in this analysis. A number of additional ‘ 
constraints were imposed on the mission, including the following: 

1. The angle between the spacecraft spin axis and the spacecraft sunline is to be between 30 
and 150 degrees for the first six months. 

2. The spacecraft is to be in sunlight during the last 15 hours of the transfer orbit. 

3. The spacecraft is to be visible from the earth at the time of retro firing. 

4. The spacecraft shadowing in lunar orbit is to be less than 2.5 hours per orbit revolution. 

Analysis showed that the first three constraints result in launch time restrictions, while the 
shadow time, which depends greatly on retro rocket firing time, trajectory dispersions, and orbit 
evolution, is not amenable to simple analysis. 


FLIGHT PLAN 


The basic elements of the flight plan are depicted in Figure 1. In the general description of 


the flight plan and vehicle given in the paragraphs below, the events and sequence are correct; 


(1) 

( 2 ) 


(3) 
(3) - ( 4) 
(5) 


( 6 ) 


(7) 

( 8 ) 

(9) 

( 10 ) 



LIFTOFF 

FIRST STAGE CUTOFF, 
SECOND STAGE 
IGNITION 

SECOND STAGE CUTOFF 
COAST PERIOD 
THIRD STAGE BURNOUT 
(INJECTION POINT) 

THIRD STAGE - SPACECRAFT 
SEPARATION 
(6 )-(7) TRACKING AND ORBIT 
DETERMINATION 
(TRANSFER ORBIT) 

RETRO IGNITION 
RETRO BURNOUT 
(INSERTION POINT) 

LUNAR ORBIT 
EMPTY STAGE HI 
CONTINUES ON 
TRANSFER ORBIT 


4 5 6 

Figure I— Flight plan of AIMP-D. 


however, the specific event times are approxi- 
mate and depend on final vehicle definition and 
trajectory shaping. 

A three- stage Improved Delta (Model 
DSV-3E) will be used as the boost vehicle. The 
first stage consists of a modified Thor booster 
with three strap-on TX-33 solid rocket motors; 
its main propulsion system is a MB3 engine 
utilizing LOX and RP-1 propellants. The 
second stage is powered by an AJ10-118 engine. 
In the first and second stages, the main engines 
are gimballed for pitch and yaw control, while 
roll control is maintained by vernier engines 
and by cold gas nozzles, respectively. The 
second stage cold gas nozzles also provide full 
attitude control during the coast period. The 
third stage is a spin- stabilized solid propellant 
X-258 motor which is mounted on a spin table 
mechanism attached to the second stage. A 
shroud on the second stage encloses the third 
stage and spacecraft during the ascent in the 
atmosphere. During first and second stage 
powered flight and during the second stage 
coast period, the booster attitude and 
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sequencing is controlled by a pre-programmed autopilot. A radio guidance system which corrects 
propulsion and autopilot deviations operates during portions of the first and second stage powered 
flight. 

Upon launch the vehicle travels vertically for four seconds, during which time a roll program 
causes the vehicle to roll from the pad azimuth to the initial flight azimuth \A t j . After four 
seconds, the programmer stops the roll program and starts the booster pitch program, which con- 
sists of four sequential pitch rates varying in magnitude and duration. The pitch program is 
designed to cause the vehicle to follow approximately a ballistic flight path (i.e., zero angle of 
attack). The three strap-on solid propellant motors attached to the booster burn out 40 seconds 
after liftoff and are dropped from the booster at 70 seconds after liftoff. Main engine cutoff 
(MECO) occurs approximately 150 seconds after liftoff. Upon MECO, the main engine thrust tails 
off and the vernier engines continue burning for approximately nine seconds. 

Four seconds after MECO the second stage separates and ignites; three seconds after ignition 
the pitch program, consisting of two rates sequentially applied to the stage, is initiated. Second 
stage engine cutoff (SECO) occurs approximately 550 seconds after liftoff. 

After SECO the spent second stage, with the third stage still attached, enters a coast period. 
During this period, the second stage autopilot programmer causes the second-third stage combi- 
nation to assume the proper attitude for third stage ignition. This maneuver is necessary because 
the third stage is spin- stabilized and then flown at constant attitude. 

At the proper time during this coast period, an automatic timer ignites the spin table rockets 
on the second stage which "spins-up" the third stage and then provides a signal to separate and 
ignite the third stage. 

At third stage burnout, the payload is on an elliptic transfer trajectory to the moon. After 
the vehicle coasts for a short time, the third stage- spacecraft separation mechanism is activated, 
and a spring imparts a separation velocity of approximately 6 feet per second. The spacecraft 
then continues on the transfer trajectory, during which various ground tracking systems obtain 
range, range rate, and other data. These data are used in a real time computer program which 
predicts and updates the transfer orbit, determining the characteristics of the possible lunar or- 
bits and related data as a function of retro rocket firing time. As the spacecraft approaches the 
moon, the latest data are used to determine an optimum retro rocket firing time. The firing com- 
mand is sent either directly or through a time-delay device. If continued tracking or telemetry 
data indicates that the retro rocket did not fire, a second signal is sent immediately. Barring 
malfunctions or excessive trajectory deviations, an acceptable orbit will be achieved after the 
retro rocket burns out. 

TRANSFER TRAJECTORY DISPERSIONS 

Transfer trajectory deviations result from variations in vehicle propulsion and guidance 
system. Deviations are introduced by the first, second and third stages; however, the largest 
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contributor by far is the spin- stabilized solid propellant third stage. While the first and second 
stages are controlled by an autopilot and guidance system, the third stage is subject to uncorrect- 
able errors in the direction and magnitude of thrust and accounts for over 90 percent of the trans- . 
fer trajectory deviations. These errors are especially pronounced if the ABL X-258 motor is 
used as the third stage. The X-258 was considered in the mission analysis; however, it is being 
replaced by the UTC FW-4 motor for the actual flight. The FW-4 motor has considerably less 
velocity magnitude error (the primary cause of trajectory deviations). Table 1 shows the approxi- 
mate three-sigma transfer trajectory injection deviations of the X-258 motor. Figure 2 shows 
two "three sigma ellipses" of transfer pericynthion location resulting from the injection errors. 

The region can completely encompass the moon; in fact, for very low energy transfer trajectories, 
it can extent past the critical Hill's surface (Reference 2) or "sphere-of-influence”, thus resulting 
in a transfer for which a lunar capture is not possible. The size of this region is highly dependent 
on the energy of the nominal transfer trajectory, increasing in length in the low energy region, 
which is the region of interest because of retro velocity limitations. These large deviations neces- 
sitate determining the conditions which optimize the probability of mission success. 

Ta ble 1 INJECTION ALTITUDE = 200 n. mi. 

Transfer Trajectory Injection: Three Sigma Deviations 

of the X-258 Third Stage Motor. MEDIUM ENERGY V. = 35590 fps 



Figure 2— Three-sigma dispersions of the transfer 
pericynthion. 


Quantity 

3 a Deviation 

Velocity 

174 ft/ sec 

Elevation Flight Path Angle 

1 .1 deg. 

Azimuth Flight Path Angle 

1.1 deg. 

Altitude 

10.5 n.mi. 


PROBABILITY STUDIES 

The primary purpose of the AIMP-D Lunar Orbit Mission Study was to determine the maxi- 
mum probability of mission success (pj compatible with the vehicle's mode of flight, spacecraft 
weight, and spacecraft requirements. The basic criterion for determining whether or not a mission 
is successful is the lifetime of the lunar orbit. Although a lifetime of at least six months is de- 
sired, a shorter lifetime would satisfy most mission requirements. The proper choice of a 
trajectory thus appears to be that which yields the highest curve of lifetime probability history, 
while favoring a point of the curve corresponding to a six-month lifetime. The basic approach in 
this analysis was to start with a simplified model in which a wide range of conditions is analyzed, 
and use the resulting data for a more complicated analysis in the smaller, more interesting region 
of conditions. 
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To determine the trajectory which maximizes P s , it was necessary to combine three studies. 
These studies were aimed at determining the following: 

1. the performance capability of the vehicle, 

2. the lunar orbit dispersions about each nominal trajectory, 

3. the orbit lifetime resulting from the lunar orbit dispersions. 

Since the studies were lengthy and the time for the combined analysis was limited, each study 
was started at about the same time and accomplished generally independently of the others. Thus, 
nominal trajectory conditions for the lunar orbit dispersion study were chosen without regard to 
vehicle capability. This approach was continued until a good understanding of the separate studies 
was available. Data from the various studies were then combined, for the region of interest, to 
obtain the final results. This paper is concerned with the second and third of the study phases 
listed above. Presented below are, first, results of approximate analytic studies on lunar orbit 
dispersions, and then results of a more sophisticated Monte Carlo analysis of these dispersions. 
These studies determined the trajectory parameters which affect mission success probability, 
their relative importance, and an understanding of the effects so that the results could be extra- 
polated and the conditions optimized. 

Probability Study Variables 

In setting up the study, one of the major considerations was the choice of a set of parameters 
to be used as input or independent variables. This set must be sufficiently complete to allow all 
of the basic effects of the problem to be analyzed, but must also be as small as possible (i.e., no 
redundancies) in order to reduce the size of the study. Many such sets may be defined: one of 
them consists of the translunar trajectory injection position, velocity, attitude, time, and the retro 
velocity, i.e., h. , p., p. , v. , y., A z , <9. , 0. , DATE, t i , andv R , where the subscript i denotes 
injection conditions.* 

These eleven parameters completely specify a translunar trajectory and retro into a lunar 
orbit (retro firing time is considered as a dependent variable). However, if this set were used 
and each parameter perturbed (one at a time) four times, the number of trajectories which must 
be studied would be 4 11 . This number is clearly unacceptable and, as will be seen, unnecessary. 
The following assumptions are made: 

1. The boost trajectory from liftoff to third stage burnout, and the vehicle attitude during 
boost, lie in a plane. 

2. The third stage centerline and the third stage burnout velocity vector are parallel. 

3. The spacecraft is not rotated in space after separation from the third stage. 


*A full list of symbols is given in Appendix A. 
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Then injection latitude (p.), longitude (p.), and azimuth ^A z j are functions only of launch 
azimuth A z and the powered flight angle ' PF . Furthermore, 


Since the assumptions depict very closely the characteristics of an actual boost trajectory, the 
set of eleven parameters can be reduced to eight parameters. 

Further reductions can be made by studying the translunar trajectory. A slightly simplified 
model which neglects the long-period terms in the lunar motion will be set up. Figure 3 depicts 

the transfer trajectory with respect to the 
earth- moon plane. 

By utilizing the geometry of Figure 3, the 
injection position vector [r 4 ] and velocity 
vector [r. ] in the X'Y'Z' coordinate system 
may be obtained from h. and V. by 

[ R i] = fij t h i + "e] > 

[K] p.JN. 

where the transformation elements are 

e i. = e (A>- ^ °tr)> 

f i. 3 f (<V 7 ?> e TR> 'Xi) ' 

and 4> 0 , v , and <9 TR are as shown in Figure 3. 

The motion of the satellite is a function of [R. J , [r. ] and the position history of the moon. 

If the lunar orbit is assumed circular, then the rotation rate and the earth-moon distance are con- 
stant and the position history of the moon after injection may be considered to be independent of 
injection time and date in the X'Y'Z' coordinate system. Using the above relationships, the tra- 
jectory after injection depends on </> 0 , v , 0 TR , y. , h t , v. , and v R . Two of these elements, 4> 0 
and <9 tr , were not compatible with inputs to the computer programs and were converted by noting 
that, for trajectories which pass within some distance of the moon, the one-to-one relationships 

A> = ^0 ( r P- Pp) = ^0 ( b - ^b) - 

®TR ” ®TR ( r p’ ^p) ” ^TR ( b ’ ^b) ’ 


Z.z 1 



Figure 3— Orientation of the trajectory plane relative 
to the Earth-Moon plane. 
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can be written if the other elements are kept constant. Further, an analysis of the matched conic 
solution showed, for a range of injection altitudes on the order of one-tenth the earth radius, that 
to a very close approximation that the transfer orbit characteristics are similar if orbit energy 
is fixed. Thus, injection velocity and altitude are related, and if a reference altitude is fixed the 
above set reduces to [v, , y., r p , p p , n, and v R ] . These six elements then replace the original 
eleven elements. They were used in the remainder of the studies. 

The launch date was deleted from the set of independent parameters which govern the transfer 
trajectory by assuming the moon to be in a circular orbit. The actual orbit has an eccentricity of 
about 0.054; thus the earth-moon distance varies by about ±11,400 nautical miles from its mean 
value. A matched conic sensitivity study revealed that this variation has a second-order effect on 
the results. 


Two-Dimensional Retro Rocket Analysis 

A two-dimensional, simplified model was used in an analysis to define, in general terms, the 
major problems associated with the lunar transfer and capture (via retro rocket) phases of the 
trajectory. While the data generated in this analysis were not exact, they were sufficiently ac- 
curate for accomplishment of the objective of the analysis; i.e., to determine the approximate 
values of the transfer trajectory parameters which would result in acceptable lunar orbits, thereby 
narrowing the range of variables to be investigated in the more detailed subsequent analyses. 

The matched conic model used in this analysis is two-dimensional, viz., the transfer orbit, 
retro rocket velocity increment (v R ) , and final lunar orbit are defined to be in the plane of the 
moon about the earth. The moon’s motion about the earth is assumed to be described by a circular 
orbit of radius (r a ) . 

The problem begins at injection into the transfer orbit, which corresponds to the Delta booster 
third stage burnout. The injection parameters given are radius (r ; ) or altitude (h.) , velocity 
(v.) , and elevation flight path angle , measured from a nonrotating geocentric frame (Figure 4). 
Given these parameters and Keplerian orbit relationships, a vectorial addition can be used to de- 
fine the velocity vector with respect to the moon. By assuming that the angle between the earth- 
vehicle line and the earth-moon line is negligible, the following relations can be derived at the 
SOI (Figure 4): 


V v/<E = [( V v/® cos y\/»Y + (V® ' V v/® sin ^v/e) 2 ] 


y v/(L 


Ggi “* n/ *) - 


The point of contact of the trajectory with the SOI , which has a radius of r s , can be defined 
by an angle y 0 . Any value of y 0 can be obtained by varying the launch time while keeping the 
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transfer trajectory fixed. The impact parameter, b, is 
given by 

b = r s sin(y v ' /c -y 0 ) . 

The above relationships specify the complete lunar 
transfer orbit. 

The one parameter remaining to be defined is the 
angle /3 between the retro rocket velocity vector and the 
earth- moon line. If the retro rocket is assumed to be 
oriented along the transfer trajectory injection velocity 
vector, then this angle is 

0 - 0 TR + y. - 90° , 

and 0 TR is a function of h. , Vi , and y r 

Given the retro rocket velocity magnitude and point 
(or time) of application, the final lunar orbit can be de- 
fined. Thus, all parameters are observed to be functions 
of the three injection parameters y 0 , V R and t R . 

The parameters V y it y 0 and V R were varied; 
for each combination the retro velocity was impulsively 
introduced at various trial times along the transfer tra- 
jectory, and the resulting lunar orbit was determined. 

If one of the trial lunar orbits satisfied the success 


Figure 4— Retro rocket analysis geometry. criterion, then that combination of variables was con- 
sidered to be acceptable. Results of this portion of the 
study are presented in Figures 5, 6, and 7 which show 
the acceptable values of the impact parameter b as a function of transfer injection velocity. For 
computational purpose^ an injection altitude of 200 nautical miles is assumed. Positive b indicates 
a leading edge approach. For the purpose of these figures, a successful orbit is defined as one 
with a pericynthion altitude greater than 1,000 nautical miles and an apocynthion radius of 25,000 
nautical miles or less. The basic quantity to note from these figures is the range b of b which 
yields successful orbits. The amount of trajectory dispersion that can be tolerated is directly re- 
lated to the width of this range. Errors in b can be related to injection errors; thus these data 
can be used to establish probability trends. 


As shown in Figure 5, for a given V R and y k the limit lines converge with increasing V. . 
There is some value of v. above which no value of b exists. No value of b exists below the mini- 
mum energy case since the transfer trajectory will not reach the sphere of influence. 


8 







The effect on the range of b of varying v R is shown in Figure 5. At the lower velocities, b 
does not change significantly with v R over a wide range in v R . The major effect of increasing v R 
is to extend the maximum value of V. at which a successful orbit can be obtained. Above 4,000 
fps the effect of v R is small. A value of 4,000 fps will be used as a basis for the discussion which 
follows since it approximates the value applicable to the actual spacecraft weights of interest. 

Figure 6 shows that at a given V. the value of b increases with increasing y. . Since the 
dispersions in b are not, to the first order, dependent on y. (References 1 and 3), an increase in 
probability can be expected with increasing y 4 . It may be shown that only 0 TR and /3 are signi- 
ficantly affected by y j . Since varying <? TR varies the nominal value of b, and this may then be 
readjusted by changing the initial phase angle, the results indicate that the effect of increasing y. 
is to align v R more nearly in direct opposition to V v/(t . Thus the trend of b resulting from a 
change in y. is the same as the trend resulting from reorientation of the retro rocket. At V i 
equal to 35,500 fps and y i equal to 0 degrees, the angle of attack at the SOI ( a S0I ) is about 36 
degrees. 



Figure 5— Effect of retro velocity on the allowable 
transfer trajectory impact parameterfor successful lunar 
orbit insertion. 


Figure 6— Effect of injection flight path angle on the 
allowable transfer trajectory impact parameter for suc- 
cessful lunar orbit insertion. 
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The nominal aiming point to maximize 
the probability of obtaining a successful orbit 
should be one which yields a value of b ap- 
proximately midway in the range of b. As 
shown in Figures 5 and 6, this aiming point in 
general corresponds to a leading edge approach. 
For high y. and low V. , b becomes extremely 
large and the location of the midpoint is not 
obvious. 

The above data basically define two- 
dimensional corridors in which transfer to a 
successful lunar orbit may be achieved. By 
combining these data with lunar impact sen- 
sitivity data (References 1 and 3), an indi- 
cation of the probabilities of successful lunar 
orbit may be obtained. The sensitivity data 
relate transfer orbit injection errors to their 
effects on the lunar impact parameter. The 
application of these injection sensitivities to 
a specific transfer corridor is illustrated in 
Figure 7. 

The transfer corridor selected to illus- 
trate this method is defined by an injection 
injection velocity, v f (fps) x io - 3 flight path angle of zero degrees and a retro 

rocket velocity increment of 4,000 fps. Exami- 

Figure 7-Effect of injection velocity on the probe- nation of the corridor presented in Figure 7 
bility of successful lunar orbit. 

again shows a generally increasing corridor 
width for a decreasing injection velocity. 

From an injection velocity of 35,450 fps to a minimum transfer velocity of 35,350 fps, the corridor 
widens to 40,000 nautical miles (the value of r s considered for this study). Thus, one expects a 
greater allowable injection-impact parameter sensitivity in the lower injection energy regions. 

The "elliptic" regions presented in Figure 7 define the impact parameter dispersion for one- 
sigma variations in vehicle performance at transfer orbit injection corresponding to the velocity 
and flight path angle errors given in Table 1. The nominal aiming points assumed for the three 
injection velocities are indicated by point circles. The ellipses may be adjusted vertically by 
simply changing the launch time (phase angle) . An indication of the probability of successful lunar 
orbit may be obtained by considering the ratio of the dispersion ellipse area within the corridor 
to the total area of the ellipse. As the injection energy is decreased, the dispersion ellipses are 
elongated. Thus, the advantage of an increasing corridor width for a decreasing injection velocity 
is somewhat negated by the increase in dispersion ellipse length. For the specific corridor pre- 
sented, the one-sigma dispersion ellipses indicate that some velocity in the 35,400 to 35,475 fps 
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region exists at which the probability of lunar orbit will reach a maximum. Lower nominal V. 
are not considered since errors greater than minus one-sigma would extend below the minimum 
energy velocity. 

In summary, the two-dimensional analysis indicated that the optimum injection velocity would 
be in the low energy region and that the probability of success will be significantly affected by the 
injection flight path angle. This latter point will be discussed in detail below. 


Monte Carlo Analysis 

Having limited the variables of concern and having determined the region of interest and an 
understanding of certain trends, the probability studies were continued using more exact and de- 
tailed techniques. This portion of the probability study was accomplished with a three-body tra- 
jectory simulation Monte Carlo program. At the outset of the study, a choice was made between 
using an approximate matched conic solution to the transfer trajectory or a numerically integrated 
trajectory solution. The former, though approximate, has the advantage of requiring less machine 
time. The choice was made to use numerical integration for the following reasons: 

1. An "exact" solution would in any case have to be programmed to check the approximate 
answers, and study time was relatively short. 

2. While the risk of using an approximate technique in a new application could be taken in a 
feasibility study, it was felt that this risk could not be afforded at the time of this study 
phase. 

Familiarity with the matched conic technique and the newer asymptotic solutions (Reference 4) 
obtained during the course of this study would in all probability have resulted in another course 
of action. 

Figure 8 is a block diagram showing the basic flow of the Monte Carlo computer program. 

The program integrates forward along the trajectory from some initial reference time; for ex- 
planatory purposes this time is defined as the second stage burnout time. Nominal initial velocity, 
position, time, and attitude components are entered. In addition, the standard deviations o-e. of 
the statistically independent vehicle parameters e. which affect the initial conditions are entered, 
along with the partial derivatives of the initial conditions to variations in these parameters. 

To obtain sample initial conditions, n random Gaussian N £ are generated and combined with 
the above data. After determining the sample initial conditions, the program integrates forward 
to third stage ignition. At this point three new n a are generated and used to compute sample 
third stage thrust and attitude angles. The program then continues integrating through third stage 
burning and along the subsequent zero thrust transfer trajectory. The distance from the satellite 
to the moon is checked, and when this distance becomes less than an input value (r s ), the program 
begins determining the lunar orbit oscillating elements which would result if the nominal fourth 
stage were fired at various points along the trajectory. The program then proceeds to a subroutine 
which selects the best firing time (discussed below). In this way the program logic parallels the 
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Figure 8 — Generalized block diagram of Monte Carlo Interplanetary Trajectory Program. 

decisions which would be made during the actual flight. Upon selecting the best firing time, the 
program generates randomly perturbed fourth stage conditions and computes the resulting sample 
burnout lunar orbit osculating elements. The program then returns, generates new random num- 
bers for the initial conditions, and repeats the total cycle until the input sample size is reached. 
The statistical distributions of apocynthion radius (r a ), pericynthion radius (r p j , eccentricity, 
inclination and the conditional (r p ) distributions for those orbits with r a < r Q (i = 1 ... 10), 
where r Q are input numbers, are then determined. 

At the time of this probability study, the only reliable method available to determine lunar 
orbit lifetime was the use of several variation-of-parameter programs which numerically inte- 
grated the lunar orbit elements. To perform this integration for six-month periods on 5 to 10 
trial firing time orbits for each of 100 samples, and then to repeat this procedure for each set of 
nominal conditions, would require an enormous amount of machine time. Thus, an alternate 
parameter-to- lifetime relationship was chosen, upon which to determine the firing time for each 
sample and upon which the probability P s was calculated. For the firing time determination, two 
techniques were chosen. The first minimized eccentricity and the second minimized apocynthion 
radius while restricting the pericynthion to the region above the lunar surface. In both cases the 
probability of mission success was estimated as 
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P s - P (r a < 25, 000 n.mi . , r p > 938 n.mi .) 

The optimum probability of obtaining a lunar orbit was investigated by parameterizing the effects 
of v. , y. , T) } r and retro velocity v R . For computational purposes, the injection conditions are 
again considered at a reference altitude of 200 nautical miles. The resulting data can be applied 
to injection at any other altitude by adjusting the velocity to maintain constant transfer orbit energy. 
December 1965 was originally chosen for the analysis, although the results are applicable for any 
month. In each lunar month, days exist with equivalent injection conditions. Prior to determining 
a series of trajectories with varying r p and p p> impact trajectories were obtained for each of the 
different injection condition combinations (V ; , y if 77). The set of trajectories will be referred to 
as the impact net. 

A set of trajectories, called a miss net, with varying r p and p p (actually, impact parameter 
b and its latitude p h were used) was obtained from an impact net trajectory. To determine a 
trajectory for a given b and p b , A z and t. were permitted to vary from those values needed 
for the impact trajectory. Thus the impact trajectory for a given miss net defines the injection 
conditions v. , y i , 77 , and a reference flight time t F , while the individual miss trajectories define 
r p and p p or their equivalent b and p\. The launch azimuth and third stage liftoff times were 
incremented so that the transfer pericynthion locations encompass a grid on the impact parameter 
(b) surface (that surface which is perpendicular to the incoming approach asymptote). 

The orbit probabilities for the trajectories comprising the miss nets were determined by 
means of the Monte Carlo Interplanetary Trajectory Program. All probabilities investigated in 
this study were for a sample size of 100. To compensate for the small sample size, several ad- 
ditional and completely independent Monte Carlo studies were performed by utilizing the ITEM 
trajectory program. In all cases the results compared to within 2 to 4 percent. The fourth stage 
was trial-fired every hour for 10 hours, the initial firing taking place when the vehicle was ap- 
proximately 20,000 nautical miles from the moon. 

Initially the firing time was chosen on the basis of minimum eccentricity. The result for a 
typical impact net trajectory is shown in Figure 9. This figure shows contours of constant P s on 
the impact parameter plane defined by varying the nominal aiming point as discussed above. As 
seen, the optimum aiming point occurs for a leading edge or retrograde transfer orbit and for b 
approximately equal to 4,000 nautical miles and p b = 45 degrees. Similar results were obtained 
for other values of V i , y. , v, and v R except at very low values of v. (35,450 fps). For V. equal 
to approximately 35,550 fps and above, the optimum b varied from about 3,000 to 4,000; and p b 
varied from zero for y i equal to zero up to 60 degrees for y i = 20 degrees. At V. = 35,450 fps, 
the optimum location of b was erratic, varying from zero to 10,000 n.mi. For the very low energy 
transfer, the region of capture and dispersion is so large that the only fact known is that the 
optimum aiming point is near the moon with either a leading or trailing edge approach. These 
trends are in good agreement with the results of the two-dimensional analysis discussed earlier. 

The effects of varying V r and 77 are indicated in Figures 10 and 11. Unlike the results of the 
approximate study it is seen that v r has a significant effect on P s for v r about 4,000 fps. In 
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Figure 9— Orbit probability versus impact parameter 
location. 


Figure 11 it is seen that probability decreases 
with increasing t>. This is explained by the 
fact that the trajectory sensitivity to velocity 
error increases as v increases. A velocity 
error has two effects: (1) it changes the true 
anomaly of the transfer orbit at r c , and (2) it 
changes the flight time and thus the moon's 
position at lunar encounter. For small 17* s 
and velocities, these two effects tend to cancel 
each other vectorially. As v increases, the 
geometry becomes increasingly three- 
dimensional and the cancellation effect 
diminishes. 

The data were reanalyzed to determine 
the effect of the retro rocket firing time deter- 
mination logic and to better define the optimum. 
The probability discussed thus far is 


INJECTION CONDITIONS: 

75=35,550 fps hj=200 n.mi. 

y ; = 0 deg . n = - 5 deg . 

PROBABILITY OF r p £938 n.mi. 



Figure 10— Orbit probability effect of retro velocity. 


PROBABILITY OF 



Figure 1 1 —Orbit probability effect of relative incli- 
nation on maximum probability. 
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P s - P (r g 5 25,000 n.mi . , r p > 938 n.mi.) , 

while the firing time was determined so as to minimize eccentricity. For this probability, the 
optimum firing time would be chosen to minimize r a while maintaining r p greater than 938 
nautical miles. Applying this criterion, typical data are shown in Figures 12 and 13. The new 
criterion yields higher probabilities; however, the trends with flight path angle and V ; are identi- 
cal within the accuracy of the results. 

For a fixed value of v and v R , the maximum value of P s on each miss net was plotted versus 
injection flight path angle and an energy- related parameter. The parameter chosen was transfer 
time to impact for the impact trajectory corresponding to each miss net. This parameter had the 
effect of expanding the scale in the low energy region. The results are shown in Figure 14. The 
maximum probability is in the range of 65 to 75 hours flight time, which if converted to injection 
velocity at 200 nautical miles (approximately 35,490 fps to 35,440 fps) is seen to be in close 

V;= 35512 fps @ h s = 200 n.mi. 
y. = 5.2 deg. 

*7 = 5 deg. 


<^qo Vr — 4000 fps 



Figure 12— Probability of successful orbit (defined as elliptic orbit with r a < 25, 000 n . mi . , r p > 938 n.mi.) 

contours on the impact parameter plane. 
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PROBABILITY OF LUNAR ORBIT 
(PERCENT) 


V ; =35584 fps@ h.=200 n.mi. 
/j =2.2 deg. 

5 deg. 

V R =4000 fps 



270° 


Figure 1 3— Probability of successful orbit contours on the impact parameter plane. 


agreement with the two-dimensional analysis results. Thus the optimum probability of obtaining 
a lunar orbit occurs for a transfer energy that results in a low enough velocity in the moon's 


•LUNAR ORBIT DEFINED AS ORBIT 
WITH PERICYNTHION RADIUS >938n. mi. 



Figure 14— Maximum probability of lunar orbit versus 
flight time to impact. 


vicinity for the retro rocket to be effective 
while keeping the transfer dispersions within 
tolerable limits. In the vicinity of optimum, 
the curve is rather flat and depends on y i , 
V R , v, etc. 

It is seen that a P s of nearly 100 percent 
can be achieved if an injection flight path angle 
of 20 degrees or greater is used. Performance 
optimization studies for the Improved Delta 
vehicle showed that, for the injection energy 
of interest, the maximum value of y i achiev- 
able when the third stage is added along the 
velocity vector (an assumption inherent in the 
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results thus far) is between 5 and 6 degrees. Figure 14 shows that this represents an increase 
of 15 to 20 percent over the zero degree y i case. This significant difference immediately rules 
out the circular parking orbit in favor of an elliptic coast orbit between second stage burnout and 
third stage ignition. 

A study revealed that varying y i had only a minor effect on the transfer trajectory sensitivi- 
ties and approach velocity but that it had a major effect on the relative angle of attack of the retro 
velocity in the moon's vicinity (Figure 4). At the sphere of influence, since y i = 0. , 

a soi ~ @ ^v/c ~ ^tr + y i ~ 90 y v /<t ' 

^ a soi 

= -2 + 1-0.25 D -1.25 deg/deg . 

Increasing y i rotates v R toward direct opposition to the approach velocity, V v/C . Monte Carlo 
runs indicated that a change in a SOI had a significant effect on P s . For the cases of interest, 
a soi was about 36 degrees for y i equal to zero. 

In order to further improve the value of y i or a S0I , the effects of pitching the third stage so 
that its angle of attack a 3 is non-zero was investigated. For this case a S0I is written as 

a SOI ~ ^TR + - 90 “ y V /C ’ 

and it is noted that 

97 i av 3 

00 ~ \ cos a 3 % 0.25 deg/deg , 

and 

d a SQI _ d^TR ^^v/C 

00. ~ dy. 00. + 1 0y i 00. 

= +0.44 deg/deg . 

The positive sign indicates that a S0I decreases and probability will increase if <9. is de- 
creased in this case. Note that for increasing probability the change in y i and 0. is negative for 
this case, while it is positive for the above case with -y. varying, with a 3 equal to zero. 

When the third stage angle of attack is increased, velocity errors due to attitude errors in- 
crease and for large values of a 3 any gain in a S0I would be lost due to the increase in av. . How- 
ever, at moderate values of a 3 this effect can improve probability as follows. The deviation in 
pericynthion radius due to attitude error is given by 
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dr p d r dV ; d r p <?y. 

dB i = dv. Je~ + dy { Je~ 


Using the impulsive approximation, we have 


V.2 n V 2 + AV 3 2 + 2V 0 AV 3 cos a 3 
3 3 



AV 3 sin a 3 

+ tan V n + AV, cosa q 
°3 


= - > 0 , 


d6>; - 


dy i 

diT. 


\ AV 3 

V. sina 3 > 


AV 3 2 + V 0, AV 3 COSa 3 

V T_ 


Combining the above equations to eliminate a 3 
or < 9 , , it can be shown that 


36 . 


= 0 


which would cancel the first order effect of 06 . 
and thus reduce the transfer trajectory dis- 
persion relative to the case with a 3 equal to 
zero. In addition, increasing a 3 reduces, though 
only slightly, the effect of ^V 3 on crV x . For this 


mission the value of 


which satisfies the 


above conditions is about 5 degrees. 


Figure 15— Monfe Carlo sample points. 


Figure 15 shows the results of a typical 
Monte Carlo run. The samples of Av t and Ay. 

are shown and the orbit lifetime for each point is indicated. The lifetime information was obtained 
by a numerical integration program. The region of infinite lifetime is indicated, and approximate 
two- sigma ellipses are shown for two values of . The trend discussed above is clearly indi- 
cated by the percentage of the ellipse area within the good region. 

Figure 16 shows the variation of P s with a 3 , for y equal to 1.66 degrees and transfer time 
equal to about 72 hours. It is seen that approximately a 5 percent increase in P s is attainable 
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with an optimum a 3 of 5 degrees. The data 
clearly indicate the effect of a 3 . This figure 
also shows the effect of oV. on P s , which is 
seen to be quite significant. These data were 
obtained for a later study, and due to a number 
of differences not discussed they are not di- 
rectly comparable to Figure 14. 

Recently, an appropriate analytic technique 
for determining lunar orbit lifetime was de- 
veloped (Reference 5) and coded. This tech- 
nique, which requires about 3 seconds of com- 
puter time per six-month lunar orbit, is based 
on the work of Lidov and Williams (References 
6 and 7) for the long-period variations, and on 
a new first order approximation to the bi- 
monthly terms due to earth perturbation. This 
technique made feasible the use of lunar orbit 
lifetime calculations to determine the retro 
firing time and to determine the value of P s 
equal to p (lifetime > six months) for a wide 
variety of cases. This technique was used to 
determine the true value of P s for the tra- 
jectories of Figure 16. The results are shown 
in Figure 17. 

As a last note on probability, the actual 
launch vehicle will utilize a UTC FW-4 third 
stage motor which has an estimated o-V 3 of 
about 0.75 the value of the X-258 aV 3 . The re- 
sulting probabilities are shown in Figures 16 
and 17. The actual launch trajectory lies on 
this curve, with a true p s value of 89 percent. 

LAUNCH WINDOW ANALYSIS 

A detailed analysis was performed to deter- 
mine the trajectory variations and associated 
performance margin required to obtain launch 
windows of up to a half hour duration. This 
study included the determination of the most ef- 
fective control parameters, correction philoso- 
phy, and the trade-offs between loss in space- 
craft weight and success probability required 


V. =35603 fps hj=169.25 n.mi. 
y. =1.66deg. o0 ■, = a>l' i =1 .06 deg. 
' V R = 3600 fps 



Figure 16— Effect of attitude on orbit probability. 


Vj =35603 fps 

r j = 1 -66 deg. 
h ; = 169.25 n.mi . 
V R =3600 fps 


P, = P (LIFETIME > 6 MONTHS) ° 9 ‘ ! 106 de 9- 



Figure 17— Effect of attitude on orbit lifetime 
probability. 
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SAMPLE SIZE = 14 



Figure 18— Probability of launch in V consecutive 
days versus launch window width. 


to obtain the necessary performance margin. 
The results indicated that, at best, a 15-minute 
launch window would result in a 20 percent loss 
in P s or a 12-pound loss in spacecraft weight. 
Neither loss could be tolerated. 

Analysis of past launches revealed that 
these losses need not occur. The countdown 
proceedings of 14 Delta launchings were re- 
viewed to determine the history and causes of 
launch delays, and whether the delay was min- 
utes, one day, or several days. The results are 
indicated in Figure 18 which shows the launch 
probability P e versus launch window for launch 
periods of one, two and three consecutive days. 
It is seen that the curves are flat for the first 
six minutes, and that in order to increase P e 
significantly with a launch window, a duration 
these curves can be explained by noting that 


of over 20 minutes is required. The shape of 
if a delay occurs it requires a finite time for correction, and invariably the countdown will 
be resumed at a point prior to the stopping point. For example, if the count is stopped be- 
tween T minus one and five minutes, the countdown will be resumed at a pre- specified point, 
say T minus 15 minutes. Thus, the likelihood of delaying the launch by only one to five 
minutes was found to be negligble. 


It is also seen from Figure 18 that the probability of launching on the first day is 56 percent 
with an essentially zero launch window. With three consecutive days, the probability is 91 percent, 
and if two launch periods are available with four days in the first period and three days in the 
second period, the probability P e is better than 99 percent. These figures are based on the sample 
proportion; applying a 90 percent confidence factor, the above number is 95 percent, which is 
quite satisfactory for mission planning. 


Thus, the launch window problem was solved by employing seven launch days divided into two 
launch periods, with a small window duration of between one and three minutes on each day. The 
trajectory profile is fixed for each day. The lines of constant injection time in Figure 11 indicate 
that for these durations the variation in P s is negligible; since the trajectory is fixed with respect 
to t £ , no performance margin is needed and thus no spacecraft weight or probability loss occurs. 


CONCLUSIONS 

In summary, the analysis determined the conditions which would maximize mission success 
probability, and showed the feasibility of obtaining a lunar orbit without a midcourse correction. 

It was found that the relative inclination v should be 5 degrees or less, that retro velocity should 
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be as high as possible, and that a relatively low energy transfer trajectory with a transfer time of 
about 72 hours should be used. In addition, spacecraft orientation had a significant effect on 
mission success probability, and in order to optimize this parameter a circular parking orbit 
could not be used. Finally, it was found that with a sufficient number of launch days, a high proba- 
bility of launch could be achieved without a significant daily launch window. 
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Appendix A 


Symbol List 

A r Azimuth, angle between the projections of true north and the velocity vector on the local 

horizontal plane 

b Impact parameter, closest approach distance to moon assuming a massless moon 

h Altitude above earth's surface 

Nj Random numbers 

n Number of vehicle errors which will perturb SECO conditions 

r Radius 

r e Earth's radius 

r c Average radius of the moon's orbit above the earth 

r t Sphere of influence radius, radius from the center of the moon defining the point at which 

the reference coordinate system is changed and lunar computations commence 

[rJ Injection position vector in earth-moon plane coordinate system 

t Time 

t f Flight time, time from injection to impact on the associated impact net trajectory 

v Velocity 

V R Retro impulsive velocity 

AV 3 Third stage impulsive velocity 

a Angle of attack 

Angle between V R and earth- moon line 
y Elevation flight path angle measured from local horizontal 

y' Elevation flight path angle measured from local vertical 

y 0 Approach phase angle at sphere of influence 

€ Vehicle error which may affect SECO conditions 

i 

ae . Standard deviation of e . 

v Relative inclination, angle between the initial transfer trajectory plane and the earth- 

moon plane 

8 Orientation angle, angle between vehicle centerline and the local horizontal 
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0 pr Powered flight angle, angle between the radius vectors from the center of the earth at 

launch and at injection measured in a non-rotating frame 

<9 tr Transfer angle, angle between the radius vectors from the center of the earth at injection 

and at lunar encounter measured in a non- rotating frame 

m Longitude 

v 0 True anomaly of moon at time of injection 

p Latitude 

P b Latitude of closest approach radius vector assuming a massless moon, measured from 

ecliptic plane 

4> 0 Phase angle 

4 Orientation angle, angle between the projections of true north and the vehicle centerline 

on the local horizontal 

ojj Argument of perigee of the lunar orbit measured from the ascending node of the lunar 

orbit on the equatorial plane 

Subscripts 

a Referring to apocynthion 

i Referring to injection conditions 

L Referring to launch conditions 

p Referring to perigee or pericynthion 

SOI Referring to sphere of influence conditions 

v/« Vehicle with respect to earth 

v/c Vehicle with respect to moon 

c/® Moon with respect to earth 

3 Third stage initial conditions 
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